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Abstract 

The  physics  and  chemistry  of  clusters  is  discussed.  First,  clusters  are 
compared  with  molecules,  on  the  one  hand,  and  then  with  solids  on  the  other. 
It  is  found  that  clusters  are  an  intermediate  state,  and  therefore  of 
special  interest.  The  Hiickel  model  is  elucidated  since  this  is  the  simplest 
of  possible  semi-empirical  methods,  and  since  it  is  readily  applied  to 
clusters.  Two  kinds  of  clusters  are  discussed  in  greater  detail:  alkali- 
metal  clusters,  because  they  are  the  obvious  candidate  for  application  of 
the  Hiickel  model,  and  Buckminster fullerene,  which  is  a  C^q  cluster  arranged 
in  a  soccer  ball  geometry. 


In  this  article  we  review  some  of  the  recent  work  reported  with 
clusters  of  atoms,  concentrating  on  the  fundamental  principles  underlying 
cluster  chemistry  and  physics.  We  begin  by  comparing  clusters  with 
molecules,  on  the  one  hand,  and  then  with  solids.  Clusters  are  an 
intermediate  state  between  individual  molecules  that  make  up  gases  and 
solids,  and  therefore  of  special  interest.  The  Huckel  model  is  elucidated 
since  this  is  the  simplest  of  possible  semi-empirical  methods,  and  it  is 
readily  applied  to  clusters.  Two  kinds  of  clusters  are  discussed  in  greater 
detail:  alkali-metal  clusters,  because  they  are  the  obvious  candidate  for 
application  of  the  Huckel  model,  and  Buckminsterfullerene,  which  is  a  C^q 
cluster  arranged  in  a  soccer  ball  geometry.  We  shall  also  consider  some  of 
the  practical  consequences  of  this  work. 

Generally,  when  we  refer  to  a  cluster,  we  mean  a  group  of  atoms, 
usually  three  to  several  hundred,  that  exist  in  a  solid  state.  An  example 
is  sodium,  which  is  a  quite  stable  solid,  but  we  do  not  usually  think  of  Na^ 
or  Na^  as  being  a  common  species.  These  would  be  considered  examples  of 
clusters. 

A  special  case  is  the  class  of  clusters  of  carbon.  Here  we  must 
distinguish  between  macromolecules  and  carbon  clusters.  Macromolecules  are 
usually  include  carbon,  hydrogen,  oxygen,  and  nitrogen,  among  others.  When 
we  speak  of  carbon  clusters,  on  the  other  hand,  we  mean  a  "molecule"  that 
consists  exclusively  of  carbon,  such  as  C^q  or  C^q.  (Even  this  can  not  be 
taken  too  literally,  since  one  reads  about  clusters  of  benzene  in  the 
literature.  However,  we  shall  not  consider  this  example  further.) 

Obviously  the  distinction  between  cluster,  molecule  and  solid  is  only  one  of 


convenience.  Nevertheless,  it  is  useful  to  discuss  clusters  separately 
because  they  have  some  of  the  properties  of  molecules  and  solids,  and  thus 
compose  a  state  of  matter  distinct  from  that  of  solids  or  gases. 

Clusters  and  Solids 

In  a  solid  crystal  there  are  a  huge  number  of  atoms,  which  can  be  taken 
as  infinite.  It  is  thus  possible  to  discuss  the  bulk  properties  of  the 
material  since  the  importance  of  the  individual  atom  or  electron  is  rather 
insignificant.  Conversely,  in  a  molecule  there  are  few  enough  atoms  that  it 
is  no  longer  possible  to  ignore  the  individual  behavior.  Let  us  look  at  the 
vibrational  motions  as  an  example.  Consider  a  non-linear  molecule  with  N 
number  of  atoms.  Each  atom  can  move  in  three  directions,  x,  y  or  z,  meaning 
that  there  are  3N  degrees  of  freedom.  In  other  words,  we  must  specify  3N 
numbers  before  we  can  describe  the  motion  of  all  the  atoms  in  the  molecule. 
However,  we  can  make  the  problem  somewhat  simpler.  First,  suppose  that  the 
whole  molecule  is  defined  to  be  at  rest.  We  can  do  this  by  using  the 
position  of  one  atom  as  the  origin  of  this  coordinate  system.  Then  by 
definition,  this  atom  cannot  move.  We  have  thus  eliminated  three  degrees  of 
freedom  called  translational  degrees  of  freedom.  Further,  suppose  we 
prevent  the  molecule  from  rotating  around  its  origin.  Since  we  have  also 
assumed  it  to  be  non-linear,  this  assumption  eliminates  an  additional  three 
degrees  of  freedom,  known  as  rotational  degrees  of  freedom.  We  are  thus 
left  with  3N-6  degrees  of  freedom.  These  are  called  vibrational  degrees  of 
freedom  in  that  they  describe  the  vibrational  motions  of  the  system. 

In  the  case  of  linear  molecules,  there  are  really  3N-5  vibrational 
degrees  of  freedom  since  there  are  only  two  rotational  degrees  of  freedom. 


The  simplest  possible  molecule  is  one  with  two  atoms.  There  are  a  total  of 
six  degrees  of  freedom  --  three  translational,  two  rotational  and  one 
vibrational  motion.  These  are  illustrated  in  Fig.  1. 

[Insert  Figure  1  about  here] 

For  molecules  it  is  relatively  easy  to  calculate  the  various  degrees  of 
freedom,  to  write  down  the  motions  of  each  atom  in  terms  of  the  so-called 
"normal  modes"  of  the  system,  and  then  to  do  explicit  calculations  of  the 
molecular  motion.  Normal  modes  are  the  fundamental  motions  of  the  system  in 
terms  of  which  all  other  possible  motions  can  be  expressed.  As  the  system 
gets  larger,  this  calculation  rapidly  becomes  impractical.  Certainly  for  a 
solid,  with  billions  of  atoms,  it  becomes  impossible  and  even  unnecessary, 
since  we  can  then  take  advantage  of  the  fact  that  we  can  ignore  individual 
atoms. 

When  confronted  with  large  systems,  it  is  no  longer  possible  to  draw 
figures  like  Fig.  1.  Instead  of  calculating  the  normal  modes,  we  now  assume 
that  the  atoms  are  arranged  in  some  periodic  repeating  structure,  called  a 
unit  cell.  We  can  then  calculate  the  motions  of  the  unit  cell,  which  may 
contain  up  to  20  atoms,  and  duplicate  that  for  the  entire  system.  But  we 
also  have  to  consider  motions  that  involve  several  unit  cells.  In  this  case 
ve  have  what  is  known  as  a  collective  motion  of  the  solid.  A  good  example 
is  a  rubber  band.  If  one  stretches  it  and  then  plucks  it,  it  will  vibrate. 
Ve  can  describe  this  motion  without  knowing  the  precise  behavior  of  each 
atom.  This  is  called  the  collective  motion  of  the  rubber  band.  For 
vibrational  motions,  which  we  have  been  discussing,  these  collective  motions 
are  known  as  phonons. 


While  we  have  been  considering  vibrational  motions,  we  could  just  as 
easily  have  been  considering  electronic  motions.  For  small  atoms  it  is 
possible  to  write  down  the  quantum  mechanical  Schrodinger  equation  for  each 
electron  and  solve  the  result  numerically.  As  the  number  of  atoms  (and 
electrons)  increases,  this  becomes  impractical,  and  so  we  instead  write  down 
the  result  in  terms  of  collective  motions.  In  the  case  of  electrons,  this 
collective  motion  is  known  as  a  plasmon. 

The  case  of  clusters  falls  in  the  middle,  however.  Clusters  are 
generally  large  enough  that  a  precise  calculation,  as  for  a  small  molecule, 
is  not  possible.  However,  they  are  still  small  relative  to  a  bulk  solid, 
and  hence  it  can  be  dangerous  to  discuss  collective  motions  in  a  cluster. 

The  study  of  clusters  is  thus  very  interesting  from  a  theoretical  standpoint 
since  they  form  a  link  between  the  macroscopic  (bulk  solid)  and  microscopic 
(molecular)  points  of  view. 

The  pictures  we  have  painted  are  at  two  opposite  extremes.  In 
actuality,  our  model  for  clusters  is  somewhere  between  an  exact  calculation 
(based  on  quantum-mechanical  first  principles  and  known  as  an  ab  initio 
calculation)  on  the  one  hand,  or  a  purely  phenomenological  calculation  (that 
treats  matter  as  if  it  were  a  continuum)  on  the  other  hand.  In  reality 
there  are  several  Intervening  possibilities  which  are  characterized  as  semi- 
empirical.  These  sorts  of  calculations  start  with  quantum  mechanical 
principles,  but  are  then  simplified  by  making  some  approximations.  These 
approximations  usually  involve  insertion  of  parameters  into  the  calculation 
that  are  chosen  solely  because  they  reproduce  the  experimental  data,  and 
hence  the  name  semi-empirical. 


Semi -empirical  methods  can  be  quite  complicated  and  also  quite 
accurate.  Their  big  advantage  is  that  they  can  be  readily  applied  to 
relatively  large  systems,  and  while  still  retaining  some  of  the  mathematical 
rigor  of  quantum  mechanics,  they  yield  relatively  accurate  results  in  a 
short  period  of  time.  One  of  the  most  popular  of  these  methods  is  known  as 
MNDO,  which  stands  for  Modified  Neglect  of  Diatomic  Overlap.*  This  method, 
available  as  a  series  of  computer  programs,  can  be  used  on  the  ubiquitous 
VAX  computer  (among  others),  and  can  easily  calculate  parameters  such  as  the 
heat  of  formation  for  molecules  containing  up  to  40  atoms,  an  impossible 
task  for  an  ab  initio  calculation. 

The  Hiickel  Model 

The  simplest  and  one  of  the  first  of  these  semi-empirical  methods  is 
2  3 

called  the  Hiickel  method,  which  was  first  developed  (in  the  1930' s)  to 
describe  conjugated  systems,  i.e.  molecules  that  can  be  written  as  an 
alternating  series  of  double  and  single  bonds  (see  Fig.  2).  In  what 
follows,  we  shall  describe  the  simple  Hiickel  method  and  then  demonstrate  its 
application  to  cluster  chemistry.  To  illustrate  the  fundamental  idea  behind 
the  Hiickel  model,  suppose  we  have  a  group  of  identical  atoms  in  a  molecule. 
For  each  atom  we  can  distinguish  between  the  valence  electrons  and  the 
inner-shell  electrons.  From  the  standpoint  of  chemical  bonds  we  can  treat 
the  inner-shell  electrons  by  assuming  that  the  valence  electrons  feel  some 
average  charge,  that  is,  the  valence  electrons  are  attracted  to  the  nucleus 
but  repelled  from  the  surrounding  electrons.  Hence  the  total  charge  felt  by 
a  valence  electron  will  be  the  difference  between  the  attraction  and  the 
repulsion,  and  this  can  be  taken  as  a  constant,  called  a,  also  known  as  the 


Coulomb  Integral.  On  the  other  hand,  a  can  be  chosen  as  zero  since  where  we 
define  the  zero-point  of  energy  is  arbitrary,  and  it  is  sensible  to  think  of 
isolated,  individual  atoms  as  having  no  energy. 

The  bond  between  the  atoms,  however,  is  due  to  the  interaction  between 
two  valence  electrons,  one  on  each  of  two  neighboring  atoms.  Since  the 
atoms  are  all  identical,  this  interaction  can  also  be  considered  a  constant, 
called  (3,  also  known  as  the  resonance  integral.  This  parameter  is  in  units 
of  energy,  and  as  long  as  one  is  not  picky  about  which  units  are  used,  it  is 
possible  to  simply  set  3  =  -1.  It  is  the  only  parameter  required  by  the 
Hiickel  model.  Since  we  are  assuming  that  the  atoms  are  bonded,  it  follows 
that  the  energy  of  this  interaction  must  be  exothermic,  and  hence  3  must  be 
a  negative  number.  The  parameter  3  tells  us  the  strength  of  the  interaction 
between  two  electrons  on  two  neighboring  atoms  and  is  thus  related  to  the 
bond  strength.  We  can  easily  relieve  our  condition  of  all  identical  atoms 
by  choosing  different  3's  for  each  species. 

[Insert  Figure  2  about  here] 

The  obvious  question  is  what  happens  when  there  are  two  valence 
electrons  on  one  atom?  In  this  case  the  Hiickel  model  no  longer  works  since 
we  have  no  way  to  account  for  this  interaction.  Electrons  on  neighboring 
atoms  are  relatively  far  away  from  each  other  and  hence  interact  relatively 
weakly.  Under  this  circumstance  it  is  possible  to  consider  an  average 
electronic  configuration  and  call  that  3*  But  if  the  electrons  are  on  the 
same  atom,  this  is  no  longer  true  and  the  Hiickel  model  breaks  down.  How 
then  can  we  apply  the  Hiickel  model  to  carbon  systems  with  conjugated  bonds, 
since  carbon  has  four  valence  electrons? 
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The  key  to  this  question  is  the  notion  of  molecular  orbitals.  Bonds  in 
conjugated  systems  consist  of  two  types,  o-bonds  and  n-bonds;  o-bonds  are 
localized  between  two  atoms,  where  each  carbon  atom  contributes  one  electron 
to  each  of  three  o-bonds,  and  the  bond  is  made  by  the  neighboring  atom 
contributing  the  second  electron,  it-bonds,  on  the  other  hand,  are 
delocalized  in  that,  in  a  fully  conjugated  system,  the  n-electrons  move 
around  the  whole  molecule.  The  essential  principle  is  that  each  carbon  atom 
contributes  only  one  n-electron,  and  furthermore  that  the  u-electrons  are  of 
significantly  higher  energy  than  the  o  electrons.  Hence  we  can  think  of  the 
11-electrons  as  being  the  valence  electrons  for  the  molecule,  whereas  the  o- 
electrons  are  thought  of  as  inner-shell.  Thus  p  stands  for  the  energy  of 
interaction  between  two  ir-electrons  on  neighboring  atoms. 

For  some  simple  systems  this  method  works  quite  well.  By  making  other 
additional  approximations,  it  is  possible  to  calculate  bond  lengths,  charge 
densities  and  other  physically  measurable  quantities.  Table  1  shows  some  of 
the  results  of  these  calculations  compared  with  experimental  data.  Given 
the  simplicity  of  the  technique,  it  is  amazing  that  the  results  are  as 
accurate  as  they  are.  It  should  be  pointed  out  that  the  actual  calculation 
is  quite  simple  and  can  be  carried  out  for  small  molecules  on  a  micro¬ 
computer.  The  procedure  is  described  in  full  in  standard  texts  on  quantum 
chemistry. 

[Insert  Table  1  about  here] 

The  Hiickel  model  has  several  implications.  First,  all  atomic 
attributes  are  expressed  in  the  constant  p.  Once  p  has  been  determined,  the 
specific  nature  of  the  atoms  is  irrelevant.  Secondly,  in  its  simplest  form, 


the  Hiickel  model  is  independent  of  bond  length.  Only  the  fact  of  whether  or 
not  atoms  are  bonded  is  of  significance.  (A  modification  of  the  simple 
technique  allows  one  to  calculate  bond  lengths.)  Therefore,  once  3  has  been 
determined,  the  results  of  the  Hiickel  model  depend  ONLY  on  the  topology  of 
the  molecule  and  on  nothing  else.  Topology  is  like  geometry,  only  simpler. 
The  topology  of  a  molecule  depends  only  on  which  atoms  are  bonded  to  each 
other.  It  makes  no  difference  how  long  each  bond  is,  what  the  bond  angle 
is,  or  how  close  or  how  far  apart  the  atoms  are.  Figure  3  contains  examples 
of  different  topologies  for  three  and  four  atom  molecules.  Recalling  that  N 
is  the  number  of  atoms  in  the  molecule,  we  can  express  the  topology  in  a 
very  simple  way  called  the  adjacency  matrix,  a  square,  N  by  N  matrix  with 
zeros  along  the  diagonal.  The  ij -component  (the  number  in  the  i1*1  row  and 
the  j1"*1  column  of  the  matrix)  is  unity  if  atom  i  is  bonded  to  atom  j; 
otherwise  it  is  zero.  An  example  of  an  adjacency  matrix  is  given  in  table 
2.  Once  we  know  the  adjacency  matrix  and  3,  then  the  Hiickel  results  are 
determined. 

(Insert  Figure  3  about  here] 

(Insert  Table  2  about  here] 

The  fact  that  the  Hiickel  model  depends  only  on  geometry  is  both  a 
blessing  and  a  curse.  Obviously  if  one  is  Interested  in  phenomena  that 
depend  sensitively  on  the  precise  configuration  around  an  individual  atom, 
or  the  precise  length  of  a  bond,  then  the  Hiickel  model  is  not  very  useful. 
If,  on  the  other  hand,  one  is  interested  in  comparing  the  effects  of 
different  geometries,  then  the  Hiickel  model  is  very  useful  indeed. 


Alkali  Metal  Clusters 

One  of  the  difficulties  of  cluster  chemistry  is  that  many  of  the 
clusters  are  very  short  lived,  and  hence  it  is  difficult  to  determine  the 
geometry  of  a  given  species.  In  many  cases,  the  only  evidence  for  the 
existence  of  such  clusters  is  mass  spectroscopy,  which  is  a  technique  which 
allows  the  accurate  measurement  of  the  molecular  weight  of  cluster  ions. 

But  it  tells  nothing  about  cluster  structure  and  geometry.  While  we  have 
considered  conjugated  carbon  systems  in  our  description  of  the  Hiickel  model, 
we  could  just  as  easily  have  considered  alkali  metal  systems,  each  of  which 
contributes  only  one  valence  electron  to  the  molecular  system.  In  fact, 
alkali  metal  clusters  have  been  an  ongoing  area  of  research  for  many  years, 
both  from  an  experimental  and  theoretical  point  of  view,  and  it  is  here  that 
the  Hiickel  model  can  be  used  with  effectiveness. 

Let  us  consider  small  alkali  metal  clusters,  ranging  in  size  from  three 
to  ten  atoms.  What  is  their  geometry?  Given  that  there  are  only  a  few 
atoms,  it  would  seem  possible  to  do  a  more  sophisticated  calculation.  The 
problem  is  that  there  are  a  huge  number  of  different  possible  geometries  to 
check.  In  principle  (and  here  we  acknowledge  that  many  of  the  possible 
structures  are  physically  unreasonable),  there  are  two  ways  of  connecting 
three  atoms,  six  ways  to  connect  four  atoms  (see  Fig.  4),  853  ways  to 
connect  7  atoms,  261,080  ways  to  connect  9  atoms  and  11,716,571  ways  to 
connect  10  atoms!  Clearly  it  is  impossible  to  do  complex  calculations  to 
determine  the  most  obvious  structure  for  all  but  the  simplest  of  molecules. 

Recently  research  has  been  done  to  calculate  the  most  stable  structure 
for  alkali  metal  clusters,  both  in  neutral  and  cationic  form,  with  up  to 


nine  atoms.  The  principle  is  to  consider  all  possible  structures  from  graph 
theory,  eliminate  those  which  were  physically  impossible  (the  computer  did 
this),  and  then  do  a  Huckel  calculation  to  determine  which  was  most  stable 
(the  computer  did  this  too).  The  results  are  shown  in  Pig.  4.  The  second 
most  stable  structures  have  also  been  calculated  but  are  not  presented  here. 
With  this  information  it  should  be  possible  to  make  some  predictions  about 
which  structures  these  clusters  assume,  and  also  to  narrow  the  field  down 
sufficiently  so  as  to  make  more  sophisticated  calculations  possible. 

Similar  research  has  also  been  carried  out  for  larger  alkali  metal  clusters, 
ranging  in  size  from  10  to  14  atoms,  albeit  not  by  going  through  every 
possible  geometry.  Instead,  the  workers  did  a  Huckel  calculation  for  likely 
geometries,  i.e.  for  geometries  that  existed  for  other  types  of  molecules. 

(Insert  Figure  4  about  here] 

The  most  interesting  results  from  this  research  are  the  conclusions 
that  the  authors  draw  concerning  the  behavior  of  these  clusters.  There  is  a 
close  similarity  between  the  Huckel  model  results  and  those  of  the  spherical 
jellium  model.  Without  going  into  the  details,  the  jellium  model  is  one  of 
the  macroscopic  ways  of  looking  at  a  solid,  which  is  very  commonly  used  in 
solid  state  physics.  However,  it  would  not  be  expected  to  apply  very  well 
to  such  small  clusters,  and  the  fact  that  it  does  is  surprising. 

Furthermore,  this  research  shows  that  there  are  fairly  distinct  breaks  in 
the  pattern.  Small  clusters,  six  atoms  or  less,  are  arranged  in  a  primarily 
two-dimensional  way.  Clusters  from  7  to  11  atoms  in  size  seem  to  have 
predominately  three-dimensional  geometry.  Clusters  larger  than  12  atoms  in 
size  contain  "body"  atoms,  that  is,  atoms  that  are  entirely  surrounded  by 


neighbors  and  do  not  lie  on  the  surface  of  the  cluster.  Hence,  one  can 
conclude  that  clusters  with  12  or  more  atoms  begin  to  behave  more  and  more 
like  bulk  solids. 


Buckminster fullerene 

Buckminsterfullerene  (Bucky  ball)  is  a  cluster  named  after  R. 

Buckminster  Fuller,  a  famous  architect  who  designed  the  Geodesic  Dome,  and 

who  incidentally  was  appointed  in  1962  as  Charles  Eliot  Norton  Professor  of 

Poetry  at  Harvard  University.  It  has  the  geometry  of  a  soccer  ball  (see 

Fig.  5)  with  a  carbon  atom  at  each  vertex.  There  are  60  identical  carbon 

atoms,  where  each  atom  is  a  part  of  one  pentagon  and  two  hexagons,  and  each 

2 

atom  is  connected  by  three  bonds  in  a  sp  hybridization.  There  are  12 
pentagons  and  20  hexagons  in  the  structure. 

[Insert  Fig.  5  about  here) 

Bucky  ball  was  first  discovered  by  researchers  at  Rice  University^ 
using  mass  spectroscopy.  When  they  first  did  the  experiment,  they  found  a 
range  of  carbon  clusters  with  sizes  from  about  32  atoms  up  to  about  80  atoms 
or  more  (the  limit  of  measurement).  But  only  even  numbers  of  carbon  atoms 
appeared  to  form  clusters;  there  were  no  clusters  with  an  odd  number  of 
atoms.  A  typical  mass  spectra  from  one  of  these  experiments  is  shown  in 
Fig.  6. 

(Insert  Fig.  6  about  here] 

The  researchers  then  allowed  the  clusters  to  react  as  much  as  possible. 
They  did  this  by  prolonging  the  time  before  the  clusters  entered  the  mass 
spectrometer,  and  also  by  mixing  free  radicals  (reactive  chemicals)  into  the 
chamber.  They  then  found  that  all  clusters  essentially  disappeared  except 


for  and,  to  a  much  lesser  extent,  C^q.  Thus  C^q  appears  to  be  a 
remarkably  stable  species.  The  question  then  arises  as  to  what  accounts  for 
this  stability. 

To  answer  that  question  we  digress  for  a  moment  and  discuss  some 
geometry.  Euler's  rule  states  that  for  any  solid  polygon  the  number  of 
faces  plus  the  number  of  vertices,  minus  the  number  of  edges,  must  equal  2. 
Stated  symbolically, 


F  +  V  -  E  =  2  . 


For  a  cube,  for  example,  F  =  6,  V  *  8  and  E  *  12.  For  a  triangular  pyramid, 
F  *  A,  V  *  5  and  E  *  7.  The  reader  can  verify  that  the  rule  holds  for  a 
square  pyramid  or  any  other  figure  (it  must  be  a  solid  polygon  without 
holes;  donut- like  shapes  do  not  follow  this  rule). 

The  most  stable  structure  for  carbon  is  that  of  benzene,  a  hexagon. 

2 

Benzene  is  stable  because  of  the  sp  hybridization  between  the  carbon  atoms, 
with  the  extra  valences  occupied  by  hydrogen  atoms.  Pure  carbon  is  most 
stable  in  the  form  of  graphite,  or  a  large  number  of  benzene  rings  in  a 
plane.  But  at  the  point  where  the  graphite  plane  is  truncated,  there  are 
loose  valences  which  would  tend  to  make  the  molecule  unstable  (unless  there 
are  hydrogen  or  other  atoms  to  occupy  the  extra  valences).  From  this 
reasoning,  it  makes  sense  to  think  that  the  pure  carbon  clusters  form  closed 
spheres  so  that,  as  in  Bucky  ball,  all  valences  are  occupied.  The  benefit 
of  this  is  that  there  are  no  loose  valences,  but  the  cost  is  an  increase  in 
strain  energy  due  to  the  fact  that  the  molecule  has  to  curve  around  on 
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itself.  Since  a  structure  of  only  hexagons  is  planar,  it  is  impossible  to 
make  a  spherical  structure  that  contains  only  hexagons. 

If  a  structure  made  up  only  of  hexagons  is  impossible,  then  the  next 
most  stable  structure  for  carbon  is  the  pentagon.  The  strain  energy  is 
somewhat  larger  in  this  case,  but  pentagonal  structures  are  known  to  exist. 
If  we  permit  a  molecule  with  both  pentagons  and  hexagons,  and  insist  that  it 
be  a  closed  polygon,  then  Euler's  rule  applies,  and  we  can  derive  the 
following: 

Let  n^  be  the  number  of  hexagons,  n^  the  number  of  pentagons,  and  N  the 
total  number  of  atoms.  Then  the  total  number  of  faces,  F,  is  the  sum  of 
pentagons  and  hexagons,  the  number  of  vertices  is  the  same  as  the  number  of 
atoms,  and  the  number  of  edges  is  “N.  The  latter  number  is  due  to  the  fact 
that  each  atom  has  three  bonds,  but  that  each  bond  is  shared  amongst  two 
atoms.  We  can  then  rewrite  Eq.  1  as 

nc  +  n,  +N-^N*nc+n,-rN*2  .  (2) 

jo  l  jo  L 
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Another  condition  is  the  limitation  on  the  total  number  of  atoms.  Each  atom 
is  a  member  of  one  pentagon  and  two  hexagons,  i.e.  three  "gons".  There  are 
five  atoms  in  a  pentagon  and  six  atoms  in  a  hexagon.  Putting  this  together 
we  get 
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Using  Eqs.  2  and  3  to  solve  for  n^  and  n^,  we  can  easily  derive  the 
following: 


n6  *  2  *  n^  *  12  .  (4) 

The  astonishing  fact  here  is  that  no  matter  how  large  N  is,  if  we  insist 
that  there  be  only  hexagons  and  pentagons  in  the  polygon,  then  there  can  be 
ONLY  12  pentagons  --  no  more  nor  less.  (The  reader  may  like  to  show  that  if 
one  constructs  a  polygon  consisting  only  of  hexagons  and  triangles,  then 
there  will  be  only  4  triangles.  The  reader  can  also  show  that  structures 
containing  only  hexagons  and  heptagons,  or  hexagons  and  octagons,  are 
impossible. ) 

The  other  interesting  fact  is  that  for  the  structure  to  be  closed, 
there  must  be  an  even  number  of  atoms.  It  is  experimentally  observed  that, 
for  clusters  larger  than  about  32  atoms,  only  even-numbered  clusters  exist! 
This  gives  strong  supporting  evidence  that  these  clusters  do,  in  fact,  form 
closed  shells.  Otherwise  there  is  no  reason  why  the  odd  numbered  clusters 
should  not  exist. 

Clearly  the  Huckel  model  should  work  very  well  for  Bucky  ball.  The 
cluster  is  highly  conjugated,  with  one  n-electron  for  every  atom.  In  fact, 
a  Huckel  calculation  of  Bucky  ball  does  show  that  it  is  energetically 
stable,  a  fact  later  substantiated  by  MNDO  calculations,^  with  further 
analysis  in  terms  of  vibrational  motions. ^  But  it  should  be  pointed  out 
that  so  far  the  only  evidence  for  the  existence  of  Bucky  ball  is  mass 
spectroscopic  (though  there  is  a  very  recent  article  describing  an  X-ray 
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crystallographic  study  that  verifies  its  existence).  No  one  has  managed  to 
synthesize  it  in  the  laboratory,  nor  isolate  it  in  a  bottle,  though  if  the 
theoretical  calculations  are  correct,  that  should  be  possible. 

Nevertheless,  the  discoverers  of  Bucky  ball  have  proposed  a  possible 
situation  of  where  it  might  occur  naturally. 

Fire  oxidizes  carbon.  The  basic  reaction,  using  methane  as  an  example, 
is 


CH4  +  202  -♦  C02  +  2H20  .  (5) 

However,  most  fires  do  not  burn  very  cleanly.  Some  of  the  material  is 
simply  vaporized  and  leaves  the  scene  as  smoke.  In  particular,  individual 
carbon  atoms  may  be  vaporized  and  eventually  form  soot.  These  individual 
atoms,  or  small  clusters  of  atoms,  are  inherently  unstable  and  reactive,  and 
so  they  would  tend  to  grow  in  size,  a  process  known  as  "upclustering. "  Many 
would  become  very  large  and  would  become  the  soot  particles  that  accumulate 
in  chimneys.  But  suppose  that  in  individual  cases  these  particles  did 
manage  to  come  together  in  a  closed  sphere.  Then,  as  Bucky  balls,  they 
would  be  relatively  stable  and  continue  up  the  chimney  as  an  element  of 
smoke. 

Conclusion 

We  have  discussed  the  current  status  of  cluster  chemistry  and  have  also 
considered  some  of  the  techniques  which  can  be  used  to  learn  more  about 
clusters.  We  have  considered  two  important  types  of  clusters:  alkali-metal 
clusters  and  carbon  clusters. 
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FIGURE  CAPTIONS 


A  diagram  illustrating  the  degrees  of  freedom  for  a  diatomic  molecule. 
The  molecule  can  move  in  any  of  three  dimensions,  x,  y  and  z.  In 
addition,  it  can  rotate  in  the  plane  of  the  page,  and  also  in  a  plane 
at  right  angles  to  the  page  (can  you  see  how  a  non-linear  molecule 
would  have  another  rotational  degree  of  freedom?).  Finally,  the 
molecule  has  one  vibrational  motion. 

Examples  of  conjugated  molecules.  Carbon  atoms  occupy  each  of  the 
vertices.  In  the  case  of  benzene,  the  hydrogen  atoms  are  shown 
explicitly.  In  the  case  of  naphthalene  they  are  not  shown  (though  they 
are  still  there),  and  the  carbon  atoms  are  numbered.  The  adjacency 
matrix  for  this  molecule  is  shown  in  Table  2. 

All  possible  graphs  that  can  be  drawn  from  groups  of  three  or  four 
atoms.  Each  of  these  structures  has  a  different  topology  since  the 
number  of  bonds  connecting  the  atoms  are  different. 

The  most  stable  alkali-metal  clusters  as  determined  by  Hiickel 
calculations.  These  clusters,  nine  atoms  and  fewer,  were  found  by 
showing  that  all  other  possible  graphs  were  either  of  higher  energy  or 
were  geometrically  Impossible.  "M"  refers  to  an  alkali-metal  atom, 
with  the  subscript  denoting  the  stoichiometry.  The  other  notation 
denotes  the  symmetry  of  the  molecule  and  is  too  complicated  to  describe 
here  (see  texts  on  inorganic  chemistry). 


A  drawing  of  Buckminsterfullerene.  There  are  carbons  at  each  vertex, 
and  no  other  atoms  in  the  molecule.  For  a  three-dimensional  model, 
look  at  a  soccer  ball. 

A  typical  mass  spectrum  for  carbon  clusters.  The  peak  height  is 
approximately  proportional  to  the  quantity  of  that  particular  cluster 
detected.  Note  that  is  the  most  prominent.  Other  experiments  have 
demonstrated  that  is,  in  fact,  much  more  stable  than  any  of  the 
others.  Note  also  that  odd-numbered  clusters  are  non-existent. 

(This  was  reproduced  from  the  first  paper  listed  in  Ref.  5  with 
permission  of  the  authors.) 
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TABLE  1 


Bond  Length  in  A 


Bond 

1-2 

2-3 

1-9 

9-10 

Modified  Huckel 

1.382 

1.414 

1.420 

1.419 

Experimental 

1.368 

1.414 

1.422 

1.419 

The  bond  lengths  of  naphthalene  as  measured  by  experiment  and  as 
calculated  by  a  modified  version  of  the  Huckel  model.  The 
numbering  system  is  shown  in  Fig.  2  (taken  from  Ref.  2,  p.  146). 
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0100000010 
1010000000 
0101000000 
0010000001 
0001000001 
0000101000 
0000010100 
0000001010 
1  0  0  0  0  0  0  1  0  1 

0  0  0  1  1  0  0  0  1  0 

TABLE  2.  The  adjacency  matrix  for  naphthalene  shown  in  Fig.  2.  Note  that 

the  matrix  is  symmetric  around  the  main  diagonal.  Renumbering  the 
molecule  will  lead  to  a  different  matrix,  but  the  Hiickel  results 
will  be  the  same. 
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